Three-loop results in HQET 
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Recent results and methods of three-loop calculations in HQET are reviewed. 



1. Off-shell HQET propagator diagrams 

Two-loop HQET propagator diagrams were re- 
duced pQ to two master integrals, using integra- 
tion by parts [2] identities. 

There are 10 generic topologies of three-loop 
HQET propagator diagrams (Fig. [I}. They can 
be reduced [3], using integration by parts rela- 
tions, to 8 master integrals (Fig. [2|). The algo- 
rithm has been constructed by hand, and imple- 
mented as a REDUCE package Grinder [3J. It is 
analogous to the massless package Mincer [4]. 

The first 5 master integrals can be easily ex- 
pressed via r functions. The next two ones are 
two-loop integrals with a single d-dependent in- 
dex. The last one is truly three-loop. 
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Figure 3. Two-loop propagator diagrams 



The first two-loop diagram (Fig. [3]) in coordi- 
nate space is just a double integral, and can be 
easily calculated for arbitrary m...s [3]: 

J(ni,n 2 ,n 3 ,n 4 ,n 5 ) = 

r(m + n 2 + 113 + 2(n 4 + n 5 - dj) 

r(m + n 3 + 2n 4 - d)T(d/2 - n i )T(d/2 - n s ) 
T(ni + n 2 + n 3 + 2n 4 — d)T (n^T (n 5 )T (ni + n 3 ) 
ni, d — 2n§, m + n 3 + 2n 4 — d 
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The second diagram in Fig. [3] for arbitrary n 
has been calculated [5] using Gegenbauer polyno- 
mials in x-space [6]: 
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Similarly to the two-loop case [7] , some HQET 
propagator integrals are related to on-shell mas- 
sive propagator integrals by inversion of Eu- 
clidean integration momenta: 
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Figure 1 . Topologies of three-loop HQET propagator diagrams 
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Figure 2. Master integrals for three- loop HQET propagator diagrams 
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In particular, at c? = 4 we have 



= -5Cs + 12C2C3 ■ (4) 

Reducing the left-hand side to the master inte- 
grals, one can express [8] the last master integral 
in Fig. [2] via the right-hand side of ((4]), which is 
an on-shell master integral known from [9|. Note 
that 0(e) corrections to ((U) are not known. 

Using this technique, the HQET heavy-quark 
propagator has been calculated up to three 
loops [10j . and the heavy-quark field anomalous 
dimension (obtained earlier by a completely dif- 
ferent method [11]) has been confirmed. The 
anomalous dimension of the HQET heavy-light 
quark current has been calculated [10]. The 
correlator of two heavy-light currents has been 
found, up to three loops, including light-quark 
mass corrections of order m and m 2 [12] . The 
quark-condensate contribution to this correlator 
has been also calculated up to three loops [12] . 
Its ultraviolet divergence yields the difference of 
twice the anomalous dimension of the heavy- 
quark current and the that of the quark con- 
densate, thus providing a completely indepen- 
dent confirmation of the result obtained in [TU] , 
The gluon-condensate contribution has been cal- 
culated up to two loops [12] (at one loop it van- 
ishes). 

2. On-shell HQET propagator diagrams 
with mass 

2.1. Reduction 

The two-loop on-shell HQET propagator inte- 
gral with a massive loop can be expressed via F 
functions for arbitrary indices [13] : 
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Figure 4. Topologies of on-shell HQET propaga- 
tor diagrams with mass 

There are two generic topologies of three-loop 
on-shell HQET propagator diagrams with a mas- 
sive loop (Fig. 2]). Algorithms of their reduction 
to master integrals, using integration by parts 
identities, have been constructed [13] by Grobner 
bases technique P3]. All scalar integrals can be 
divided into two classes which are not mixed by 
recurrence relations: apparently even and appar- 
ently odd ones. They are integrals which would 
be even or odd with respect to v —v if there 
were no iO in denominators. 

All apparently even integrals of the first topol- 
ogy reduce to 




while apparently odd ones to 

^ ^ 

All apparently even integrals of the second topol- 
ogy reduce to 




A.G. Grozin 



while apparently odd ones to 



The master integrals 




can be easily expressed via T functions. The mas- 
ter integral 




has been investigated in detail |15|16j . 

2.2. A master integral 

Now we shall discuss the integrals 
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(/in is one of the master integrals). Several ap- 
proaches have been tried |13|17j . The best result 
was obtained [T7] using a method similar to [15] . 
First we consider the one-loop subdiagram 



(7) 
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If po < 0, we can deform the integration contour 
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This integral is 



7rV2r(n2) 

(fc 2 - l)(d-l)/2-na 



dfc- 



(2fc - 2 Po y 



Inm 2 (po) = 

r(ni + n 2 - 2 + e)r(ni + 2n 2 - 4 + 2e) 
r(n 2 )r(2(m+n 2 -2 + e)) 



/: . , ni,ni + 2n 2 - 4 + 2e 
2 1 1 ni+n 2 - | +e 



or, after using a 2 Fi identity, 
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This result was obtained p~7] using the HQET 
Feynman parametrization. 

Now we can integrate in d d ~ 1 p in the three-loop 
diagram: 
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The square of 2 Fj in (|8|) can be expressed via 
an 3 F 2 using the Clausen identity. We analyt- 
ically continue this 3 F 2 from 1 + p EQ > 1 to 
z = 1/(1 + p E0 ) < 1 and integrate © term by 
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term. The result contains, in general, three 4F3 
of unit argument. 

The convergent integral /122 is related to the 
master integral Jm by 



(rf-3) 2 (rf-4)(3d-8)(3d- 10) 
122 ~ 8(3d- ll)(3d- 13) 

For this integral, we obtain [T7] 
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Expansion of this result up to e 7 agrees with [T7] 
I122 7r 2 r 3 (l + 2e)r 2 (l + 3 £ ) 
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This equality has also been checked by high pre- 
cision numerical calculations at some finite e val- 
ues. During the workshop, David Broadhurst has 
rewritten this conjectured hypergeometric iden- 
tity in a nice form 
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We have no analytical proof. 



2.3. Other master integrals 

Now we turn to other master integrals [13] . 
By applying Mellin-Barnes techniques to the a- 
representation and using the Barnes lemmas, we 
succeeded in calculating one integral exactly: 



T(l/2 - e)r(-e)r 2 (2e)r(l + e)r(3e - 1) 

4r(3/2-e)r(4e) 
[t/;(1/2 - e) + V(l - e) - 2 log2 + 2 7E ] . (13) 

In other master integral calculations, we use 
the Mellin-Barnes of the massive two-point func- 
tion: 
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Applying this representation and eliminating 
one of the lines in the triangle with unit indices 
using integration by parts, we obtain a single 
Mellin-Barnes integral: 
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Figure 5. Topologies of three-loop massive on-shell propagator diagrams 
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Similarly, in this case, we apply (| 14[) once and 
5C ©: 
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Feynman integrals considered here were 
used [13] for calculating the matching coeffi- 
cients for the HQET heavy-quark field and the 
heavy-light quark current between the 6-quark 
HQET with dynamic c-quark loops and without 
such loops (the later theory is the low-energy 
approximation for the former one at scales below 
m c ). Another recent application — the effect 
of m c 7^ on b — > c plus lepton pair at three 
loops [T5]. The method of regions was used; the 
purely soft region (loop momenta ~ m c ) gives 



integrals of this type. Two extra terms of e ex- 
pansion of the master integral of Sect. 12.21 were 
required for this calculation which were not ob- 
tained in [13] . This was the initial motivation 
for PI]. 

3. On-shell massive QCD propagator dia- 
grams 

These diagrams are used for calculation of 
QCD/HQET matching coefficients P2]. There 
are two generic topologies of two-loop on-shell 
propagator diagrams. They were reduced [201151 
121] to three master integrals, using integration by 
parts. 

There are 11 generic topologies of three- loop 
HQET propagator diagrams (Fig.JSJ) (10 of them 
are the same as in HQET, Fig. [TJ plus the dia- 
gram with a closed heavy-quark loop which does 
not exist in HQET). They can be reduced [11] 
to 19 master integrals (Fig. [6]) using integration 
by parts. The reduction algorithm was first im- 
plemented as a FORM package SHELL3 [TT]; in 
some papers [22], reduction is performed using 
the Laporta method. The master integrals are 
known from [9|llj (see also [22]). 

As an example, let's consider chromomagnetic 
interaction of the heavy quark in HQET. The 
two-loop anomalous dimension of the chromo- 
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Figure 6. Master integrals for three-loop massive on-shcll propagator diagrams 
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magnetic interaction operator has been calculated 
in [23J, using off-shell HQET diagrams (Sect. [J) 
and R* operation to get rid of infrared diver- 
gences. It was confirmed |24j . a week later, using 
a completely different approach — QCD/HQET 
matching. In the same paper, also the coeffi- 
cient of the chromomagnetic interaction opera- 
tor in the HQET Lagrangian has been calculated 
up to two loops. Recently, both the anomalous 
dimension and the interaction coefficient have 
been calculated at three loops [25], also using 
the QCD/HQET matching. These results con- 
tain two colour factors which don't reduce to Cf 
and Ca (due to gluonic "light-by-light" diagrams; 
in the /3-function and basic anomalous dimensions 
of QCD, these colour factors first appear at four 
loops. 
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